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Studies of threshold production of top quark pairs will 
provide information about Standard Model parameters 
and a window to "New Physics" phenomena. Because 
the t quark decays fast compared to the strong interac- 
tion time scale, even the threshold top pair production 
in e + e~ or 77 collisions can be described using perturba- 
tive QCD PJ. However, close to the threshold, the cross 
section is enhanced by the strong attraction between the 
quark and the antiquark and an all-order resummation 
of the Coulomb effects is required. The magnitude of 
the threshold cross section and the position of its peak 
are sensitive to the top quark mass and width, as well as 
the strong coupling constant a s , and can help determine 
those parameters with high precision 0,0. I n addition, 
precision measurements in the threshold region at the 
Next Linear Collider (NLC) will probe electroweak inter- 
actions of the t quarks and be sensitive to "New Physics" . 

A particularly interesting environment for studying 
"New Physics" effects in the top production is the 77 
collider mode of the NLC, with laser beams backscat- 
tered from the high-energy e + and e~ as the source of 
the colliding photons j| . It has been pointed out || |(| 
that the process 77 — ► tt is twice as sensitive to anoma- 
lous photon-top quark couplings as e + e~ — > tt, and in 
addition is free from contamination by the Zti coupling. 
It will therefore allow searches for anomalous couplings 
of the top quark to photons, in particular the electric 
dipole moment. Although precision studies might be fea- 
sible only at the NLC, the observation of 77 — * tt may 
be possible already at the LHC M. 

Although the threshold production cross section can be 
computed perturbatively, it is not clear how well the se- 
ries converges. Recently, the renormalization group tech- 
nique was employed || to sum up large logarithmic cor- 
rections appearing in threshold problems. It was applied 



to a(e + e~ — > tt) at the threshold in j^] and an achieve- 
ment of a 3% accuracy in the normalization of the cross 
section was claimed. However, it is conceivable that pos- 
sible large non-leading-logarithmic effects may spoil this 
picture. To a large extent, study of the tt threshold pro- 
duction cross section in the 77 mode might help in an as- 
sessment of the situation and provide an additional test 
of the assumptions of the resummation program. 

In general, behavior of the tt threshold cross sections 
is rather similar in e + e~ and 77 collisions; the difference 
comes from the fact that the quarks are mainly produced 
in spin one and spin zero state in e + e" and 77 collisions, 
respectively. Since the Coulomb interaction does not de- 
pend on the spin, this fact does not affect the behavior 
of the leading order cross section but, since relativistic 
corrections do depend on the total spin of the produced 
fermion pair, it becomes important in higher orders. 

There is yet another interesting feature of the top 
quark threshold production in 77 collisions In con- 
trast to an e + e~ collider, manipulating polarization of 
the incoming photons can suppress the 5-wave produc- 
tion and provide a possibility to study the P-wave thresh- 
old production without the huge S'-wave background. All 
these features of the top quark threshold in 77 collisions 
make it a very interesting laboratory for studying the 
non-relativistic QCD dynamics of a heavy quark anti- 
quark system, determining the standard model parame- 
ters, and searching for "New Physics" . 

A known problem with this program is the monochro- 
maticity of the photon beams at the NLC where a typical 
energy spread is currently estimated to be about 10% of 
the total energy (see e.g. jH^). This may hamper the 
threshold studies since a large energy spread will wash 
out any pronounced signal at the top threshold. It re- 
mains to be seen if this difficulty can be overcome. 
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Much effort has already been spent on theoretical 
studies of 17(77 — * ti) with polarized photons. Away 
from the threshold, one-loop QCD corrections were com- 
puted in jlM. Close to the threshold, the cross section 
<r(e + e~ — > tt) has been studied in NLO in B. Recently, 
studies of NNLO corrections to cr(e + e~ — > tt) were com- 
pleted and the appropriate effective field theory frame- 
work for studying higher order corrections to the thresh- 
old phenomena was formulated (see |3j and references 
therein) . Applying similar techniques to 77 — » tt, it is in 
principle straightforward to obtain the NNLO corrections 
to this process; the major obstacle has been the lack of 
the two-loop matching coefficient for an effective opera- 
tor responsible for the leading order transition 77 — * tt 
at the threshold. In this paper we describe the calcula- 
tion of this matching coefficient and present the NNLO 
corrections to 17(77 — > tt) at the threshold. 

Close to the threshold, the produced top quarks are 
non-relativistic so that it is reasonable to expand the pro- 
duction amplitudes in powers of the relative velocity (3 
of t and i. The Born amplitude for 77 — > ti with 0(f3 2 ) 
accuracy can be written as 



M 



-8inaQt(i 



[M s + M P ]x, 



(1) 



where Qt = 2/3 is the electric charge of the top quark, 
and x are the top and anti-top spinors. The amplitudes 
that give rise to the production of the ti state in S and 
P waves are, respectively 




(jm) 2 jf 



in ■ lei x eM 



to 2 2m 2 , 
(jm)(an)(eie 2 ) + (pei)(ae 2 ) + {pe 2 )(aei) 



(2) 



In these formulas, n is a unit vector along the photon 
flight direction, p is the top quark three-momentum and 
ei 2 a re the polarization vectors of the colliding photons. 
From now on we will consider the incoming photons to 
be circularly polarized. If the two photons have opposite 
polarizations, the S wave production amplitude vanishes. 
As a consequence, only M. p contributes to the top pro- 
duction and the cross section at the threshold behaves as 

(7-| ~ /3 3 . If the photon helicities are the same, the top 

quarks are mainly produced in the S wave with a small 
admixture of the P wave appearing in the NNLO. 

It is customary to define a normalized unpolarized 



cross section, R~ 



(7(77 



ti)/ao, with <7q 



Ana /(3s), and decompose it into contributions with the 
same and opposite photon helicities, i? 77 = (R ++ + 
i? H )/2, where we have used R ++ = R and i?" 1 = 
R valid for electromagnetic processes. Let us con- 
sider the case of the two photons with equal helicities 
and select the contribution of the S wave. (The case 
of the P wave production has been investigated earlier 
in the literature (for the most recent discussion see fl2| ) 
and, because of its relative 0{f3 2 ) suppression, we do not 
have much to add to this issue.) In case of the S wave, 



the quarks are produced in a spin singlet state. The ex- 
pression for the cross section in the regime a s <C (3 <C 1 
reads: 



Rl 



= 6QtN c (3 



,(3) 



where j3 — y/l — 4m 2 / 's. Also, 



A M 



7T~ 



+ Y P + 0((3 2 ). (4) 



The two-loop corrections can be decomposed into four 
parts, arising from the abelian and non-abelian gluon ef- 
fects, and to light quark and top quark vacuum polar- 
ization insertions in the one-loop correction. Neglecting 
terms of O {(3), the two-loop corrections are 

A( 2 ) = C F A A + C A A NA + T R N L A L + T R N H A H , (5) 
5 I 



A„ = 



12(3 2 
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+ f 7r2 + f + ||- 4 -2- 2 ln(2/?) + 2^; 



+7T A I - - H 2ft) + 2x NA ; 
A L = (-^- + lln(2p))^- + 2x L 



18 



P 



2x H . 



(6) 



In the above formulas the terms x a ,xn a ,xl,xh are 
related to the hard renormalization of the operator 
<f> + n ■ [e*i x e*2] X responsible for 77 — * ti transition at 
the threshold: 



x A = -21.02; 
xl = —0.565; 



xna = -4.79; 
x H = 0.224. 



(7) 



To evaluate the hard renormalization factors we 
have used some of the results obtained for the para- 
positronium decay (l3| [l4|. In addition, to find xna 
which is our main new result here, we had to evaluate 
the non-abelian diagrams (see Fig. Ill(a-e)) and massless 
insertions in the one-loop diagrams (an example is shown 
in Fig. 0(f)). Here we briefly summarize those calcula- 
tions. Diagram ^(a) is the simplest two-loop non-abelian 
contribution. It is both ultraviolet (UV) and infrared 
(IR) finite and we evaluate it using Monte Carlo inte- 
gration over six Feynman parameters. (We use FORM 
p"5| for symbolic manipulations and Vegas (l6| for nu- 
merical integrations.) Diagram 0(b) is UV divergent. To 
compute it we use a trick: assign mass M to the quark 
line (inside the loops) and expand the diagram in the 
ratio t = m 2 /M 2 , treating M as much larger than m. 
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FIG. 1: Examples of non-abelian hard corrections to 77 — > 
computed in this paper. In addition, we include diagrams 
similar to (f) with gluon, Faddeev-Popov ghost and fermion 
loops inserted in all one-loop corrections. 



At the end we have to evaluate the sum of the resulting 
series at t = 1. It is sufficient to compute only a few 
first terms of the series in t if we change the variables, 
t — > — 4w/(l + w) 2 , expand in w, and evaluate the new 
series at w = 2\/2 — 3 ~ —0.17. 

Diagrams (c) and (d) are also UV-divergent, with a 
subdivergence in the non-abelian vertex correction and 
an overall divergence. It is useful to first remove the 
subdivergence by subtracting a product of one-loop di- 
agrams which can be evaluated analytically. Next, the 
remaining overall divergence is removed by subtracting 
a similar two-loop amplitude which depends either only 
on the external quark momentum (for (c)) or only on the 
spatial photon momentum (for (d)). The resulting finite 
integrals are evaluated with Vegas, while the subtracted 
divergent pieces can be computed either exactly (for (c)) 
or with the expansion used for the diagram (b). 

The only IR-divergent non-trivial two-loop diagram 
is (e). It also contains a UV-subdivergence which we 
remove in the same way as in (c) and (d). The IR- 
divergence is subtracted by neglecting the loop momen- 
tum in the "t-channel" quark propagator. That subtrac- 
tion can be computed analytically, since it is simply a 
2-loop vertex diagram at threshold, of the type for which 
we developed a general algorithm in an earlier study [fl7| . 
The finite difference is evaluated with Vegas. 

We did not include effects of the 77 — > gg scattering 
via light quark boxes. We have checked with a rough 
approximation that this finite and gauge- invariant subset 
contributes only insignificantly. 

Using the result (||) it is easy to obtain the top thresh- 
old production cross section for the photons with equal 
helicities close to the threshold. In doing that, we follow 
the approach described in Ref. |Q. Let us first present 
a general formula for the resummed cross section for two 



R s + 



2inN c Qf 



different choices of the photon helicities. We find: 

1 

-C h 

\In^ 11-^/3* ) G sins (0,0, E + iT t )) , (8) 



5 ,2 



where Ch is the hard renormalization factor and G slng 
is the Green's function of the spin singlet state of the ti 
pair including relativistic corrections and corrections to 
the Coulomb potential (see |ll| for details). 

As we already mentioned, there is also a P wave contri- 
bution to the ++ cross section at NNLO. It is suppressed 
by O(0 2 ) relative to the 5 wave and working to NNLO 
we need to know it only to leading order. Since for the 
photon helicities H — the top pair is also produced in the 
P wave, there is a leading-order relation between i?p + 
and i? H , 



R+- 
R+- 



- 3 R P , 
32nQfN c 



(9) 



<9 x <9 y Im G (x, y, E + iT t ) | x=0 , y =o- 



Numerically (see e.g. [[12]) the P wave contribution is 
small; the corresponding values of i? H are ~ 5 x 10~ 2 . 
For this reason, we do not consider the P wave contribu- 
tion in what follows. In Fig. || we present the LO, NLO 
and NNLO excitation curves for Rg + computed using the 
pole mass of the top quark. As is clearly seen from this 
plot, the NNLO corrections to the normalization of the 
cross section are quite large: close to the peak they are 
about 20 — 30%. The position of the peak of the cross 
section from which the mass of the top quark is to be 
determined suffers from significant shifts when one goes 
from LO to NNLO. All these features are quite similar 
to the known behavior of e + e~ — > ti at the threshold [[|. 

Let us now address these problems in turn. It was ar- 
gued Jl9], |2Cj in connection with the threshold production 
of ti in e + e~ collisions that significant shifts in the po- 
sition of the peak are the consequences of the fact that 
the pole mass scheme is unstable against radiative correc- 
tions. The way out of this problem is to adopt a different 
mass definition which will have such a stability. In Fig. |3| 
we show the 5-wave part of the cross section, parame- 
terized by the so-called kinetic mass We observe 
that the stability of the peak does improve significantly, 
a behavior familiar from e + e~ — > ti studies. 

Let us now discuss the corrections to the peak height. 
Again, the situation is similar to e + e~ — > ti but the cor- 
rections are larger. Since the height of the cross section 
is controlled by the wave function at the origin, ip(0), of 
the ti ground state, it is easy to understand what is hap- 
pening by looking at the lna s enhanced corrections to 
ip(0) computed in [Q. Neglecting effects of the running 
of the coupling constant, they read 

ip 2 {0) = ipl(0) [l-a^lna s (7.55- 1.19 5(5 + 1)) 

+a 3 s In 2 a s (-9.48 + 0.99 5(5 + 1))] . (10) 
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FIG. 2: R+ + at LO, NLO and NNLO. Top quark pole mass 
m t = 175.05 GeV width Y t = 1.43 GeV are used. The hard 
renormalization scale and the factorization scale are equal to 
the top quark mass (a s (m t ) = 0.109). The soft renormaliza- 
tion scale is 60 GeV (a s (60 GeV) = 0.127) but we also show 
the NNLO curve for the soft renormalization scale 30 GeV to 
demonstrate large changes with the coupling constant scale 
variation. 
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FIG. 3: R-s + parameterized with the kinetic mass 
?7ikin(15 GeV) = 173.10 GeV. Other parameters are the same 
as in Fig. ^. Dashed, dashed-dotted and solid curves are LO, 
NLO and NNLO approximations, respectively. 

The 0(al) (N 3 LO) leading logarithmic effects are in- 
cluded in this formula. From Eq. (|l^) one sees that there 
is, effectively, a compensation between spin-dependent 
and spin-independent terms in the correction to the wave 
function at the origin; since the top quarks are produced 
in the spin-singlet state in 77 — > tt, the corrections in 
this case are expected to be larger than in e + e~ — > tt. 
Note also that in both cases the a 3 In a s corrections are 
negative, so that the NNLO curves shown in Fig. || are 
expected to be pushed down in N 3 LO, which clearly il- 
lustrates the sign-alternating nature of the perturbative 
series. The success of the resummation program [^| for 
e + e~ — > tt is related to this property of the series; for this 
reason one can expect a significant improvement in the 
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FIG. 4: The ratio of Rg + and R e in the threshold region at 
LO, NLO and NNLO (with two different soft renormalization 
scales) parameterized by kinetic mass mki n (15 GeV). 



stability of the height of the cross section also in 77 — * tt 
once the resummation program is carried out. 

It is also interesting to consider the ratio of the top 
threshold production cross sections in e + e _ and 77 col- 
lisions, plotted in Fig. [|. We see that the ratio is energy 
independent at LO and NLO, since in these orders the 
spin of the tt pair decouples. However, at the NNLO the 
non-relativistic Hamiltonian depends on the spin of the 
tt pair and the ratio of the cross-sections becomes energy 
dependent. This results in large corrections to the ratio 
in the vicinity of the peak. Away from the peak where 
the bound state dynamics is not very important and the 
ratio is, essentially, given by the hard renormalization 
factors of the production currents, the convergence of 
the perturbative series for the ratio is quite good. 

The NNLO corrections to the top pair threshold pro- 
duction in 77 collisions, which we have completed by 
evaluating the non-abelian hard renormalization factor, 
turn out to be large. This is similar to the process 
e + e~ — > tt. In both cases we see that the position of 
the cross section peak is stabilized when a short-distance 
quark mass definition is used, but the height of the peak 
is changed significantly by the NNLO corrections. We 
have argued that the structure of the NNLO corrections 
to 77 — * tt at the threshold is rather similar to that of 
e + e~ — ► tt: in particular the normalization of the cross 
section is determined by sign-alternating series. This 
suggests that if the renormalization group improvement, 
applied to e + e~ — > tt in || is applied to 77 — > tt, one 
can expect significant improvement in the stability of the 
normalization of the production cross section. For this 
reason, in the future, it would be interesting to resum 
the logarithmic corrections O (ln(/3)) for 77 — > ti, along 
the lines of such an analysis for the e + e~ annihilation 
H . We hope that a combined analysis of both processes 
will shed new light on the behavior of higher-order cor- 
rections to the threshold processes and allow to reach 
better theoretical precision. 
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